Color-magnitude relations of late-type galaxies 



Ruixiang Chang 1 ; Shiyin Shen l ; Jinliang Hou 1 ; Chenggang Shu 2,1 ; Zhengyi Shao 



O 
O 

(N 



O 

(N 



> 

in 
o 

O 

a: 

6 



ABSTRACT 

We use a large sample of galaxies drawn from the Sloan Digital Sky Survey (SDSS) and 
Two Micro All Sky Survey (2MASS) to present Color-Magnitude Relations (CMRs) for late- 
type galaxies in both optical and optical-infrared bands. A sample from SDSS Data Release 4 
(DR4) is selected to investigate the optical properties. Optical-infrared colors are estimated from 
a position matched sample of DR4 and 2MASS, in which the photometric aperture mismatch 
between these two surveys is carefully corrected. It is shown that, after correcting the dust 
attenuation, the optical colors for faint galaxies (i.e. M r > —21) have very weak correlation 
with the luminosity, while the optical colors for bright galaxies (i.e. M r < —21) are redder 
for more luminous galaxies. All (optical, optical-infrared and infrared) colors show similar but 
stronger correlations with stellar mass than with absolute magnitude. The optical colors correlate 
more strongly with stellar mass surface density than with stellar mass, while optical-infrared and 
infrared colors show stronger correlations with stellar mass. By comparing the observed colors of 
our sample galaxies with the colors predicted by stellar population synthesis model, we find that 
massive late-type galaxies have older and higher metallicity stellar population than less massive 
galaxies. This suggests that CMRs for late-type galaxies are trends defined by the combination of 
stellar mean age and metallicity. Moreover, our results suggest that the stellar mean metallicity 
of late-type galaxy is mainly determined by its stellar mass, while the star formation history is 
mainly regulated by the stellar mass surface density. 

Subject headings: Galaxies:evolution - Galaxies:photometry - Galaxies: stellar content 



1. Introduction 

It has long been known that early-type galaxies 
in nearby clusters exhibit a tight Color-Magnitude 
Relation (CMR), i.e. more luminous early-type 
galaxies have redder colors than less luminous ones 
(Baum 1959; Visvanathan & Sandage 1977; Bower 
et al. 1992), and this finding has been confirmed 
by later studies (Aragon-Salamanca et al. 1993; 
Stanford et al. 1995, 1998; Ellis et al. 1997; Ko- 
dama et al. 1998; Terlevich et al. 2001; Bernardi 
et al. 2003, 2005; Blanton et al. 2003a; Bell et 
al. 2004; Holden et al. 2004; Hogg et al. 2004; 
Cool et al. 2006; Chang et al. 2006). This rela- 
tion has been mainly ascribed to the metallicity 
effect, either because more massive galaxies have 
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stronger binding energies so that the gas can be en- 
riched to higher metallicity (Faber 1977; Dressier 
1984; Arimoto & Yoshii 1987; Kodama & Ari- 
moto 1997), or because larger ellipticals form at 
low redshift through major mergers between more 
massive, metal-rich disk systems (Kauffmann & 
Chariot 1998; De Lucia et al. 2004; Kang et al. 
2005). 

However, the case for late- type galaxies is much 
more complicated partially because of the influ- 
ence of on-going star formation and the effect of 
dust attenuation. Almost thirty years ago, Vis- 
vanathan & Griersmith (1977) find that early-type 
spiral galaxies (SO/a to Sab) in the Virgo cluster 
exhibit very similar optical CMR to that of E/S0 
galaxies, i.e. brighter spirals show redder colors, 
except that the former has a much larger scatter. 
Based on a large sample from Sloan Digital Sky 
Survey (SDSS, York et al. 2000), Blanton et al. 
(2003b) investigate the broadband optical prop- 
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erties for concentrated galaxies (fitted Sersic in- 
dex n > 3) and exponential galaxies (n < 1.5) 
and find that exponential galaxies have a less 
tight, but undeniable optical CMR. Baldry et al. 
(2004) analyze the bimodal color-magnitude dis- 
tribution of SDSS galaxies and find that both red 
and blue-distribution galaxies can be well fitted 
by a straight line plus a tank function. In other 
words, the low-luminosity blue-distribution galax- 
ies (M r > —19) show a shallower CMR slope than 
high-luminosity galaxies. However, these works 
are all limited to the optical bands. 

Tully et al. (1998) investigate global extinction 
in spiral galaxies and establish optical-IR OMRs 
for spirals after correction for dust attenuation. 
Peltier & de Grijs (1998) determine an optical- 
infrared CMR (I — K vs. Mk) for spiral galaxies. 
They find the CMR is very tight and its slope is 
steeper than that for elliptical and SO galaxies. 
This result is consistent with early studies (e.g. 
Wyse 1982; Tully et al. 1982; Mobasher et al. 
1986). Using stellar population synthesis (SPS) 
models, they conclude that faint spiral galaxies 
have both a younger age and a lower metallicity 
than bright spirals. Bell & de Jong (2000) have 
compiled a sample of 121 low- inclination spiral 
galaxies with radially resolved optical and near- 
infrared photometry. By comparing the avail- 
able photometry with the SPS model colors, they 
found that the star formation history of a galaxy 
(as probed by its age, metallicity and gas frac- 
tion) strongly correlates with its surface bright- 
ness and magnitude. These results are confirmed 
by Macarthur et al. (2004). However, all those 
early studies are limited to small galaxy samples 
since it is difficult to get large and homogeneous 
optical and infrared photometric samples simulta- 
neously. 

On the other hand, the interpretation of these 
observed CMRs is still a challenge for semi- 
analytic models, van den Bosch (2002) and Bell et 
al. (2003) find that there is significant discrepancy 
between the predictions of semi-analytic models 
and the data. The observations show that brighter 
spiral galaxies are redder than fainter ones, while 
the models predict that fainter spirals should be 
slightly redder than bright spirals. Recently, Kang 
ct al. (2005) introduce a semi-analytic model 
based on high-resolution N-body simulations and 
compare their model predictions with the observed 



bimodal distribution in the color-magnitude dia- 
gram found by Baldry et al. (2004). They find 
that the prediction for the red branch is very 
close to the observations, but for the blue part, 
the predicted colors of luminous spirals are bluer 
and those of less luminous spirals are redder than 
observations. All these results suggest that the 
implied physics of gas accretion and feedback dur- 
ing the formation and evolution of spiral galaxies 
are still not well understood (van den Bosch 2000; 
Bell et al. 2003). 

It should be pointed out that the comparison 
presented by Kang et al.(2005) is only limited to 
the color [u — r) , which is very sensitive to dust at- 
tenuation. Furthermore, the degeneracy between 
stellar age and metallicity for optical colors com- 
plicates further investigations for the star forma- 
tion history of late-type galaxies. Thus, there is 
still a long way to go before getting robust conclu- 
sions. Since the combination of optical and near- 
infrared broad-band colors can partially break age- 
metallicity degeneracy, it is important to estab- 
lish observationally the CMRs in both optical and 
near-infrared bands for late- type galaxies based on 
a large homogenous sample. 

The SDSS and Two Micron All Sky Survey 
(2MASS, Jarrett 2000a, b) have already created 
the largest well-defined sample of galaxies with 
well-measured optical and infrared photometry to 
date. Recently, we have analyzed the optical and 
infrared color properties for early-type galaxies 
after correcting the aperture mismatch between 
these two magnitude systems (Chang et al. 2006). 
This paper is a complementary study of Chang 
et al. (2006). Here the samples are selected for 
late-type galaxies. One of our aims is to estab- 
lish optical, optical-infrared and infrared CMRs 
for late-type galaxies by using the largest homo- 
geneous samples up-to-date. 

To derive the optical-infrared colors, we adopt 
the method presented in Chang et al. (2006) and 
correct SDSS magnitudes to the aperture where 
2MASS magnitudes are measured. For the pur- 
pose of minimizing the different dust attenuation 
caused by different inclinations, the galaxies in our 
sample are limited to face-on spirals. We also cal- 
culate the emission-line-free colors by comparing 
the magnitude difference before and after remov- 
ing the emission lines from the spectra. With these 
carefully corrected optical and optical-infrared col- 
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ors, wc compare the correlations between colors 
and absolute magnitude, stellar mass and stellar 
mass surface density and try to distinguish which 
parameter mainly determines the color properties 
for late-type galaxies. Moreover, by comparing the 
colors of our sample galaxies with those of SPS 
models, we investigate the stellar populations of 
spiral galaxies and their implications on star for- 
mation processes in late-type galaxies. 

The outline of this paper is as follows. Section 
2 describes the way of selecting sample galaxies 
and the photometric quantities used in the paper. 
Section 3.1 presents the method to do the photo- 
metric aperture corrections. The magnitude cor- 
rections for emission-line contamination is shown 
in Section 3.2. Section 4 is the main part of the pa- 
per, where we present various correlations between 
colors and absolute magnitudes, stellar mass and 
stellar mass surface density and their comparisons. 
Section 5 describes the SPS models and the model 
predictions for stellar populations of spiral galax- 
ies. Finally, in Section 6 we summarize the main 
results of this paper. 

2. The data and samples 

The data used in this paper is mainly taken 
from Data Release 4 of the SDSS (York et al. 2000; 
Stoughton et al. 2002, hereafter EDR; Adelman- 
McCarthy et al. 2006). In particular, we use 
the Main Galaxy Sample in the New York Univer- 
sity Value-Added Galaxy Catalog (NYU-VAGC), 
which is based on the SDSS database, but with an 
independent set of enhanced reductions (Blanton 
et al. 2005). From this catalog, we take spec- 
troscopic redshift and photometric quantities, in- 
cluding Petrosian magnitudes, K-corrections, axis 
ratio (b/a) from the exponential fitting (a and 
b are the semimajor and semiminor axes) and 
fracDev r . The parameter fracDev r is obtained 
by taking the best-fit exponential and de Vau- 
couleurs fit to the surface brightness profile, find- 
ing the linear combination of the two that best fits 
the image, and storing the fraction contributed by 
the de Vaucouleurs fit. We also take the quantities 
from ProfMean, the azimuthally averaged surface 
brightness in a series of fixed circular annuli (sec 
table 7 in EDR) , which are used to make the aper- 
ture correction for the optical-infrared colors (see 
section 3.1). 



The 2MASS is a ground-based, near-infrared 
imaging survey of the whole sky and its Extended 
Source Catalog (XSC) contains nearly 1.6 million 
galaxies (Jarrett et al. 2000a, b). The infrared 
magnitudes are taken from extended source cata- 
log (XSC) of 2MASS. In section 3.1, wc will de- 
scribe how to estimate the optical-infrared colors. 

Conversion from apparent magnitude m to ab- 
solute magnitude M depends on the adopted cos- 
mology and the K-correction: 

M = m-5log[D L (z,n M ,n A )] -25-K(z), (1) 

where Dl is the luminosity distance in unit of 
Mpc and K(z) is the K-correction. We assume 
the standard cosmology with Am = 0.3, ft a = 0.7 
and H n =70 km s _1 Mpc -1 to calculate luminosity 
distance. We use SDSS r— band Petrosian mag- 
nitude and calculate M r at redshift z = 0. K- 
corrections for magnitudes in the g, r, i, z bands 
(and J, H, K bands if needed) are taken from 
NYU-VAGC (Blanton et al. 2003a; Blanton et al. 
2005). Stellar masses of galaxies are taken from 
http://www.mpa-garching.mpg.de/SDSS (see de- 
tails in Kauffmann et al. 2003a). 

Sample galaxies are selected from Main Galaxy 
Sample in the NYU-VAGC (Blanton et al. 2005). 
We adopt the r— band morphology indicator 
fracDev r to select late-type galaxies with cri- 
teria fracDev r < 0.5, which means that the sur- 
face brightness of selected galaxy is more likely 
to be expressed by exponential profile. Adopting 
fracDev parameter has the advantage of seeing 
and inclination correction during the profile fit- 
ting (Bcrnardi et al. 2005). From this criteria, 
our sample contains 154,870 late-type galaxies. 

We further restrict our sample with the condi- 
tion b/a > 0.75 so as to select nearly face-on galax- 
ies based on two considerations. Firstly, edge-on 
galaxies are more dust-attenuated and reddened 
than face-on systems (Tully et al. 1998; Xiao et 
al. 2006). This effect will introduce extra disper- 
sion into the CMRs. In other words, by select- 
ing face-on galaxies, we can minimize the effects 
of dust attenuation caused by different inclination 
of the spiral disks. Secondly, the isophoto aper- 
ture of the high axis-ratio galaxies is nearly cir- 
cular. In this case, the SDSS ProfMean measured 
in circular radii can be properly used to make the 
aperture correction of optical-infrared colors with 
enough accuracy. After this, the number of late- 
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apparent magnitude m r redshift z 

Fig. 1. — Distributions of several quantities for 
the galaxies in Sample A (solid lines) and Sample 
B (dotted lines), including R^{arcsec), concen- 
tration C = i?9o /R50 , r-band apparent magnitude 
m r and redshift z. 

type galaxies is reduced to 45,033. 

Furthermore, we restrict our sample to galaxies 
in the redshift range 0.01 < z < 0.2. The lower 
limit is to reduce the uncertainty using the redshift 
as distance estimator. The upper limit is to lower 
the uncertainties on K-corrections for the magni- 
tudes and to reduce the possible color evolution 
with redshift. Finally, the number of galaxies is 
further reduced to 40,987. We denote this sample 
as Sample A and use it to investigate the optical 
CMRs. 

The NYU-VAGC also contains a matched sam- 
ple between SDSS and 2MASS XSC. This sample 
is obtained by positionally matching sources ob- 
served by 2MASS to galaxies in the 'main' spec- 
troscopic sample of the SDSS within iarcsec. We 
select all the galaxies in sample A with 2MASS 
match to create Sample B. There is no cut in nei- 
ther rriK nor surface brightness in K-band. The fi- 
nal number of galaxies in Sample B is 9,930. Sam- 
ple B is used to investigate the optical-infrared and 
infrared CMRs. 

Fig. 1 plots histograms of several quantities for 
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AR/R 50 flux ratio in r-band 

Fig. 2. — Histograms of several quantities rel- 
evant to do aperture corrections for galaxies in 
Sample B, including two times Petrosian radius 
in the r-band 2Rp (from SDSS), the isophotal ra- 
dius R k2 afc (from 2MASS), AR/R 50 and the r- 
band flux within Rk20fc divided by the Petrosian 
flux. 

the galaxies in Sample A (solid lines) and Sample 
B (dotted lines), including R^{arcsec), concen- 
tration C = i?go/i?5o (-Rgo and -R50 are the radii 
including 50% and 90% of the Petrosian fluxes), 
r-band apparent magnitude m r and redshift z. As 
we can see from the top-right panel of Fig. 1, 
most of the sample galaxies have C < 2.6, con- 
sistent with the early studies of the morphology 
separation with parameter C (e.g. Strateva et al. 
2001). Compared with Sample A, the galaxies in 
Sample B are biased to be apparently bright, large 
and to lower redshifts. This is due to the fact that 
2MASS is much more shallow than SDSS. 

3. Corrections for magnitudes 
3.1. Aperture corrections 

In order to estimate the optical-infrared colors, 
it is important that both SDSS and 2MASS mag- 
nitudes are measured within the same aperture. 
We adopt the same method as that in Chang et 
al. (2006) to do aperture corrections for the mag- 
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Fig. 3. — Histograms of emission-line correc- 
tions for magnitudes in g— band (left-top), r— band 
(right-top), i— band (left-meddle) and colors g — r 
(right-middle) ,r — i (left-bottom) and r — z (right- 
bottom) in the fiber. 

nitudes, i.e. to correct SDSS magnitudes to the 
aperture where the 2MASS magnitudes are mea- 
sured. 

There are several sets of magnitudes provided 
in the 2MASS XSC. In this paper, we use the 
isophotal fiducial magnitudes, which are measured 
within the circular aperture corresponding to a 
surface brightness of 20.0 mag arcsec -2 in the K s - 
band (the aperture is denoted as Rk20{c)- SDSS 
provides the azimuthally averaged surface bright- 
ness in a series of circular annuli, i.e. ProfMean. 
We calculate g,r,i,z magnitudes that are matched 
to the 2MASS measurements by interpolating the 
cumulative radial surface brightness profile in each 
SDSS band at the corresponding isophotal radius 
Rx20fc (using the recommended approach in sec- 
tion 4.4.6.2 of EDR). The derived magnitude in r- 
band m r K is used to estimate the optical-infrared 
colors throughout this paper. The average value 
of magnitude correction in the r-band for aperture 
mismatch is 0.26. 

Fig. 2 plots the distribution of several quanti- 
ties for galaxies in Sample B, including two times 



Petrosian radius in the r-band 2Rp (from SDSS) 3 , 
the isophotal radius Rx20fc (from 2MASS), 
the aperture mismatch Ai?/i? 50 , r = (2Rp — 
RK20fc) / Rbo,r and the ratio between the r-band 
flux within Rk20}c and the Petrosian flux. It 
can be seen that the aperture mismatch between 
SDSS and 2MASS is quite large. However, more 
than half of the Petrosian flux is already contained 
within Rx20fc for most of galaxies. 

3.2. Emission-line corrections 

Since late-type galaxies have on-going star for- 
mation, emission lines from HII regions have sig- 
nificant contributions to the observed magnitudes. 
However, SPS models generally do not include 
contributions from emission lines. Therefore, it 
is necessary to correct the emission-line contami- 
nations for the magnitudes when we compare the 
predicted colors from stellar population models 
with the observed ones (as what we will do in Sec- 
tion 5). 

We use the magnitudes computed from the 
galaxy spectra to correct emission-line contamina- 
tions for the magnitudes. The magnitude differ- 
ence estimated from the spectra before and after 
removing the emission lines (denoted as AM em ) 
represents the correction inside the fiber. Assum- 
ing the contribution of emission lines to the to- 
tal magnitude equals to that in the fiber 4 , the 
sum of Petrosian magnitude and AM em can be re- 
garded as emission-free magnitude and the colors 
estimated from this set of magnitudes are treated 
as emission-line-free colors. 

Fig. 3 presents distributions of emission-line 
corrections for magnitudes and colors in the fiber. 
It can be seen that, although the corrections for 
the magnitudes are small, emission lines do influ- 
ence the colors. 



3 The r-band 2Rp is the aperture to which the Petrosian 
magnitudes are integrated, see EDR 

4 This assumption may not always be correct. Since the outer 
part of spirals has bluer colors, which implies a higher aver- 
age specific star formation rate (SFR/M*) than that in the 
inner part, our method may underestimate the emission- 
line contamination. On the other hand, due to the contri- 
bution of AGN in the nuclear region of massive galaxies, 
our method may overestimate the emission-line contam- 
ination. 
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Fig. 4. — Histograms of r — i (after emission-line 
and dust corrections) for galaxies in Sample A at 
different bins. The solid curves are the best-fit 
results for Gaussian functions. The median value 
of M r in each bin is shown as the number in each 
panel. 

3.3. dust corrections 

Although dust attenuation in face-on late-type 
galaxy is not expected to be large, it may have 
significant effect on the colors. It would be help- 
ful if the influence of dust can be estimated and 
accounted for. 

Kauffmann et al. (2003a) has already ob- 
tained the dust attenuation in the z-band, A z , 
for SDSS galaxies. They use the 4000 A break 
strength and the Balmcr absorption line index 
HSa to constrain the best star formation history 
for one galaxy. The dust extinction is estimated 
by comparing observed imaging g — r and r — i 
colors to model-predicted colors. We adopt A z 
from http://www.mpa-garching.mpg.de/SDSS for 
our sample galaxies. The median value of A z for 
Sample A and B is 0.21 mag and 0.24 mag, re- 
spectively. By adopting the standard attenuation 
curve of the form t\ oc A -0 7 , we can estimate 
the dust attenuation in other bands. We refer 
"colors" as these emission-line- free and dust- free 
colors throughout this paper if there is no special 
notation. 

It should be pointed out that our method to 
correct the effect of dust attenuation may not be 
accurate for individual galaxy. However, the re- 
sults for the global trends are robust. Indeed, 



Obric et al. (2006) find evidences to confirm that 
the estimates of A z are related to the galaxy dust 
content. 

4. Results 

4.1. Color-magnitude relations 

In this section, we use Sample A (only from 
SDSS) to investigate optical CMRs and Sample 
B (a matched sample of SDSS and 2MASS) for 
optical-infrared and infrared CMRs. The opti- 
cal colors (g — r, r — i, r — z) are estimated 
from the Petrosian magnitudes. The optical- 
infrared colors (r — J, r — K) are derived using the 
aperture-corrected magnitudes in the r-band and 
the isophotal fiducial magnitudes from 2MASS. 
The infrared color J — K are estimated from the 
isophotal fiducial magnitudes from 2MASS. All 
the magnitudes are corrected for Galactic fore- 
ground extinction and K-correction. 

The galaxies in each sample are divided into 
15 bins according to their absolute magnitudes in 
the r-band, M r . We choose equal intervals of M r 
in each bin and make sure that each bin contains 
more than 100 galaxies. The colors are corrected 
for dust extinction by calculating a median value 
of A z in each bin since the estimated A z for indi- 
vidual galaxy is model-dependent but the general 
trend is reliable. It is found that the median A z 
in luminous bin is larger than less luminous bin 
and the median A z is slightly negative for faint 
galaxies (M r > -19). 

Fig. 4 plots histograms of r — i (after emission- 
line and dust corrections) for galaxies in Sample 
A at different bins. The median value of M r in 
each bin is shown as the number in each panel. It 
is clear that the distributions of r — i can be well 
approximated by Gaussian functions, which arc 
plotted by solid curves obtained from maximum 
likelihood estimations. The distributions of other 
colors show similar behaviors. Therefore, we can 
assume that the color distributions in each bin can 
be expressed by a Gaussian function 
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which is characterized by the mean C{M r ) and 
the dispersion a(M r ). The best estimations for 
C(M r ) and a(M r ) are listed in Table 1. Further 
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Table 1: color-magnitude relations 
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tests show that main results do not change if a 
narrower redshift range (such as 0.01 < z < 0.1) 
is used to construct the sample. 

The results for CMRs are plotted in Fig. 5, 
where the filled square represents the value of the 
color (y-axis) as a function of M r (x-axis) in each 
bin, the error bar represents the dispersion a(M r ). 
We also show the results after emission-line correc- 
tions but without dust corrections as open squares 
without error bars. To quantify the tightness of 
the relation, we define a parameter S = — , where 
AI is the difference between the mean color of the 
brightest bin and the faintest bin, a is the average 
value for the color dispersions in 15 bins. Thus, 
higher S corresponds to tighter relation. The re- 
sulted S for the CMR is shown in the upper-left 
corner of each panel in Fig. 5. 

It can be seen that dust attenuation contributes 
significantly to the observed CMRs. This is ex- 
pected since luminous galaxy is more attenuated 
and reddened than less luminous object (Tully et 
al. 1998; Masters et al. 2003). After correcting 
the dust attenuation, the optical colors {g — r, r—i, 
r — z) for faint galaxies (i.e. M r > — 21) have very 
weak correlation with the luminosity, while the op- 
tical colors for bright galaxies (i.e. M r < —21) are 
redder for more luminous galaxies. This is con- 



sistent with the results of Baldry et al. (2004), 
who also find a shallow CMR slope for the low- 
luminosity blue-distribution galaxies. 

It is also shown that, after correcting dust at- 
tenuation, the optical-infrared (r—J, r—K) CMRs 
are still very tight. For infrared color (J — K), a 
weak CMR is presented, in the sense that the color 
range AI is comparable with the mean dispersion 
er. Unfortunately, due to the fact that 2MASS is 
much more shallow than SDSS, there is no suffi- 
cient number of less luminous galaxies in Sample 
B to give the optical-infrared and infrared CMRs 
behavior at the fainter end. 

4.2. Correlations between colors and stel- 
lar mass 

Stellar mass is another fundamental physical 
parameter for galaxies. Now, we turn to inves- 
tigate the color-mass relations for late-type galax- 
ies. 

Again, we divide the galaxies into 15 bins ac- 
cording to their stellar mass M* (which is in units 
of Mq throughout this paper) and fit the distribu- 
tion of the color in each bin by a Gaussian func- 
tion. The method to correct dust attenuation is 
the same as that in section 4.1. The best-fit mean 
of color and its dispersion in each bin are listed in 
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18 -20 



Fig. 5. — Color-magnitude relations, where the filled square represents the value for the emission-line-free 
and dust-free color (y-axis) as a function of M r (x-axis) at each bin. The error bar represents the dispersion 
a(M r ). Open squares without error bars are the results after emission-line corrections but without dust 
corrections. The x-axis of open square is slightly shifted for clearness. The number in the upper-left corner 
of each panel shows the result of S = — for each correlation (see text for details). 



Table 2 and plotted in Fig. 6. The most remark- 
able feature of Fig. 6 is the tightness of correlation 
between optical-infrared color (r — J, r — K) and 
stellar mass, in the sense that AJ is nearly (if not 
more than) more than 3 times the mean color dis- 
persion. We will discuss the implication of these 
results in Section 5. 

It could be interesting to compare the CMRs 
with color-mass relations. The colors show similar 
but stronger correlations with stellar mass than lu- 
minosity. The colors in the reddest bins are redder 
and the colors in the bluest bins are bluer when 
mass is used instead of absolute magnitude. This 
is true for all wavelength bands. However, it's not 



completely clear how much of the steeper slopes 
and tighter scatter in color-mass relations rela- 
tive to CMRs is driven by the color dependence 
of stellar M/L values. Tremonti ct al. (2004) also 
find that, when luminosity is used instead of mass, 
both dust and M/L variations act to smear out 
the turnover seen in the mass-metallicity relation 
at high metallicity. 

4.3. Correlations between colors and mass 
surface density 

Stellar mass surface density, which describes 
how luminous mass is distributed along the disk, 
is another fundamental parameter for spiral galax- 
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Fig. 6. — Correlations between colors and stellar mass. The caption is the same as that of Fig. 5. 



ies. The half-light radius in the z-band and stellar 
mass yield the effective stellar mass surface density 
yu* = M*/27ri?5 (Kauffmann et al. 2005; Chang 
ct al. 2006), which is in units of M kpc~ 2 . The 
results for color-/z* relations are given in Table 3 
and Fig. 7. 

It can be seen that, both optical and optical- 
infrared colors show correlations with stellar mass 
surface density, i.e. galaxies with higher /i* have 
redder colors. This point is anyhow expected since 
high mass surface density corresponds to higher 
stellar mass for late- type galaxies (Shen et al. 
2003). More interestingly, when ^u* is used in 
place of M„, the resulted S increases for optical 
color, while S decreases for optical-infrared and in- 
frared colors. In other words, for late-type galax- 
ies, the optical colors correlate more strongly with 
mass surface density than with stellar mass, while 



optical-infrared and infrared colors show stronger 
correlations with stellar mass. 

This is obviously different from what is ob- 
tained for early-type galaxies. Chang et al. (2006) 
have shown that, for elliptical galaxies, all galaxy 
colors correlate more strongly with stellar mass 
than with mass surface density, which suggests 
that the star formation history of nearby elliptical 
galaxies is primarily determined by the mass of 
the system. Since the age-metallicity degeneracy 
makes it difficult to directly deduce the property of 
stellar population from colors, it is instructive to 
break the age-metallicity degeneracy firstly, such 
as adopting the SPS models, and then discuss the 
implications of various correlations on the star for- 
mation histories (SFHs) for late-type galaxies. 
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Table 2: Correlations between colors and stellar mass. 
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Table 3: Correlations between colors and stellar surface mass density log(/i*) 
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Fig. 7.- 
Fig. 5. 



Correlations between colors and stellar mass surface density. The caption is the same as that of 



5. Stellar populations of late-type galaxies 

Motivated by the work of Bell & de Jong (2000) , 
we use a combination of optical and near-infrared 
broad-band colors for galaxies in our sample B to 
probe the stellar populations of late type galaxies. 
In other words, we generate a grid of model spec- 
tral energy distributions using SPS models and es- 
timate the model parameters by matching all the 
available photometry of the sample galaxies to the 
colors predicted by the SPS models. By doing this, 
we are able to investigate the trends in the mean 
stellar age and metallicity as a function of stellar 
mass and mass surface density. 

We adopt the SPS models of Bruzual & Char- 
lot (2003). The universal initial mass function 
is taken from Chabricr (2003). The lower and 



upper mass cut-offs are set as m^=0.1M Q and 
m{7=100M Q , respectively. We use simplified SFHs 
and given metallicities to construct a grid of model 
colors. The star formation rate (SFR) is parame- 
terized by: 



*(*) 



if t < tf orm 

Bte -t/r SF if tfQrm < t <t g 



(3) 



where B is a normalization constant which deter- 
mines the total mass of the stellar population, tsf 
is the star formation timescale, t g is the age of 
galaxy and is fixed as t g = 12Gyr, tf orm charac- 
terizes the time delay for star formation. Under 
the above assumptions, the mean age of the stellar 
population < age > is given by: 



< age >= — 



(4) 
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Therefore, the mean age of stellar population is 
determined by the combination of tsf and tf orm - 
If tform = 0, the youngest mean age given by 
none- negative t$f is around 4Gyr, which corre- 
sponds to t — ► oo and the SFR increases with 
time. In this case, the stellar population is still 
too old to fit the colors of some observed blue 
galaxies. As an alternative, we assume a none- 
zero tf orm , which insures that our models also in- 
clude the youngest stellar populations. We assume 
tform — 0.2rgF during model calculations. This 
assumption is acceptable for two reasons. First, 
we focus only on the mean stellar age rather than 
details of SFH for each individual galaxy. Fixing 
the relation between tsf and tf or m will reduce the 
number of free parameters in the model. Second, 
SPS model has the advantage of giving robust rela- 
tive ages and metallicities estimations (see discus- 
sions in Bell & de Jong 2000). We will emphasize 
the relative trends rather than the absolute values 
of ages and metallicities. Indeed, test calculations 
show that the main conclusions are not sensitive to 
either the adopted form of SFH or the assumption 

tform — 0.2r S F- 

Thus, there are two free parameters in the 
model: metallicity and tsf- t~sf progresses from 
to 60Gyr (tsf — 60Gyr corresponds to the upper 
limit for i/ orm = \2Gyr) smoothly. This results 
in the mean age of stellar population going from 
old to young. Firstly, we generate a model grid 
of single-burst stellar population (SSPs) for differ- 
ent metallicity by interpolating along the model 
metallicities. The range of the model metallici- 
ties is from 0.005Z Q to 2.5Z Q . Then, we gen- 
erate a model grid for a number of t$f values 
by integrating the SSPs with a given metallic- 
ity. These models provide a set of model colors 
Cmodei,i(< age >,Z) (where i denotes different 
color) for a range of SFHs and metallicity. 

We match all the available colors of the galaxies 
to the colors predicted by the SPS models and 
estimate the model parameters using % 2 -fitting. 
For each tsf and metallicity Z, the corresponding 
X 2 is given by: 

2 \ ^ [Cmodel.i C bs.i) / r \ 

X =2., —2 — ■ (5) 

i=i 1 

where <7j is the observed error for color i. The 
best model match (< age >,Z) is the one with 
the minimum x 2 value. 



We use the observed g — r, r — i,r — z,r — J 
and r — K colors in x 2 -fitting. These colors are 
corrected by emission-line contamination and dust 
attenuation. Since the estimated A z is model- 
dependent and may not be accurate for individ- 
ual galaxy, we choose to correct an average trend 
of dust attenuation by fitting a linear relation be- 
tween A z and M*, i.e. to estimate A z from stellar 
mass. In general, massive galaxies have larger dust 
attenuation. The results for % 2 -fitting show that 
there are about 5 percent of galaxies with the min- 
imum x 2 value located in either boundary of the 
model parameters tsf and Z. We exclude these 
galaxies from further analysis below. 

Fig. 8 plots the model predictions of (< age > 
, Z) against stellar mass and mass surface density 
for galaxies in Sample B. To show the distribu- 
tion more clearly, we plot sets of contours to indi- 
cate the positions within which different fraction 
of total number of galaxies is contained. We also 
adopt the Spearman rank-order test to quantify 
the significance of the relation and show the re- 
sulted correlation coefficient in the right corner of 
each panel. Fig. 8 shows that massive galaxies 
have older and higher metallicity stellar popula- 
tions than less massive galaxies. In other words, 
not only the stellar metallicity but also the stellar 
mean age that contribute to the observed CMRs 
for late- type galaxies. To contrast, the CMRs 
of early-type galaxies are mainly ascribed to the 
metallicity effect: more luminous galaxies have 
redder colors because they have higher metallici- 
ties (Faber 1977; Dressier 1984; Arimoto & Yoshii 
1987; Kodama & Arimoto 1997; Kauffmann & 
Chariot 1998; De Lucia et al. 2004; Kang et al. 
2005). 

The general trend for mass-metallicity relation 
found in this paper is similar to some earlier works 
(such as Peltier & de Grijs 1998; Bell & de Jong 
2000). Moreover, gas-phase oxygen abundance in 
star-forming galaxy has already been found to be 
related to the luminosity or/and stellar mass of 
the system (Kobulnicky & Zaritsky 1999; Kob- 
ulnivky & Kewley 2004; Hammer et al. 2005; 
Liang et al. 2004; 2006 etc.). For comparison, we 
show the observed mass-metallicity relation pre- 
sented by Trcmonti et al.(2004) with three blue 
lines. To transfer 12 + log{0/H) to log(Z/Z Q ), 
we assume the metallicity Z can be represented by 
oxygen abundance and take the solar value as 8.8 
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Fig. 8. — Model predictions of mean age (< age >) and metallicity (Z) against stellar mass and mass surface 
density for galaxies in Sample B. The set of contours indicates the positions within which how much fraction 
of total number of galaxies is contained. The levels of contours correspond to 68%, 90% and 95%, respectively. 
The number in the right corner of each panel shows the Spearman rank-order correlation coefficient. Three 
blue lines in bottom-left panel plot the observed rang for the relation between stellar mass and gas-phase 
metallicity, which is taken from Table 3 in Tremonti et al.(2004). The solar 12 + log(0/H) is taken as 8.8 
(Grevesse & Sauval 1998). 



(Grevesse & Sauval 1998). Since the metallicity 
derived from SPS models is the mean value of the 
stellar population, it is reasonable that the model 
predicted mettallicty is somewhat lower than that 
of Tremonti et al. (2004) . As we can see, the gen- 
eral trend of our results deduced from the OMRs 
are consistent with their independent measure- 
ments from emission lines. 

More interestingly, it can be seen (even by eye) 
that the metallicity correlates more strongly with 
M* than /i*. We also find that the residuals from 
the age-/x* relation do not correlate with M*, while 
the residuals from the age-M* relation increase 
with /«* . All these suggest that the stellar mass of 



late-type galaxy is the primary parameter to de- 
termine the mean metallicity of the system. This 
is probably, similar as early-type galaxy, because 
the stellar mass strongly correlates to the poten- 
tial well of the galaxy and then determines the 
efficiency for ejecting the metals of the system. 
Tremonti et al. (2004) also interpret their results 
as strong evidence of both the ubiquity of galactic 
winds and their effectiveness in removing metals 
from galaxy potential wells. 

On the other hand, the resulted mean stellar 
age < age > is more sensitive to /x* than . It is 
also found that the residuals from the metallicity- 
/x* relation strongly correlate with M*, while the 
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residuals from the metallicity- M* relation have no 
correlation with yu*. In other words, our results 
suggest that mass surface density plays an more 
important role in regulating the star formation his- 
tory of late-type galaxy than stellar mass. Indeed, 
some early studies already pointed out that the 
star formation history for late-type galaxies corre- 
lates more strongly with mass surface density than 
stellar mass (Bell & de Jong 2000; Kauffmann et 
al. 2003b; 2006). Kennicutt (1998) also shows 
that it is local (e.g. surface density) rather than 
global factors (e.g. stellar mass) that regulate the 
rate at which spiral galaxies form their stars at the 
present day. 

6. Summary 

We have collected two samples of face-on late- 
type galaxies to investigate the correlations be- 
tween colors and absolute magnitude, stellar mass 
and stellar mass surface density. Sample A is se- 
lected from DR4 of SDSS and Sample B is a po- 
sition matched sample between SDSS (DR4) and 
2MASS. In order to correct the aperture mismatch 
between SDSS and 2MASS, we use the radial pro- 
file of SDSS galaxies and correct the SDSS mag- 
nitudes to the isophotal circular radius where the 
2MASS magnitudes are measured. The emission- 
line contamination is corrected for the magnitudes 
by comparing the magnitudes computed from the 
spectra before and after removing the emission 
lines. We also correct dust attenuation by adopt- 
ing A z estimation plus an standard kind of attenu- 
ation curve. Moreover, we match the observed col- 
ors of our sample to the colors predicted by stellar 
population synthesis model and investigate their 
implications of the stellar populations. We find 
that: 

1. The optical colors for faint galaxies (i.e. 
M r > — 21) have very weak correlation with 
the luminosity, while the optical colors for 
bright galaxies (i.e. M r < —21) are red- 
der for more luminous galaxies. The optical- 
infrared colors show tighter CMRs than op- 
tical ones. 

2. All (optical, optical-infrared and infrared) 
colors show similar but stronger correlations 
with stellar mass than with absolute magni- 
tudes. 



3. The optical colors correlate more strongly 
with mass surface density than with stellar 
mass, while optical-infrared and infrared col- 
ors show stronger correlations with stellar 
mass. 

4. By comparing the observed colors with SPS 
model colors, we find that massive galax- 
ies have older and higher mctallicity stellar 
populations than less massive galaxies. This 
suggests that the CMRs of late- type galaxies 
are the combined results of stellar mean age 
and metallicity. 

5. It is also shown that the mean stellar metal- 
licity is more sensitive to than to 
while the mean stellar age is more sensitive 
to [i* than M*. This suggests that the stellar 
metallicity is mainly determined by the mass 
of the system, while the star formation his- 
tory of late-type galaxy is mainly regulated 
by mass surface density. 

Our results are based on the homogeneous and 
large galaxies samples from SDSS and 2MASS to 
derive the CMRs for late-type galaxies in both op- 
tical and optical- infrared bands. Since the com- 
bination of optical and near-infrared broad-band 
colors contains more information on the star for- 
mation process, our observed results have provided 
strong constraints on the model descriptions for 
gas accretion and supernova feedback processes in 
late-type galaxies. Further comparisons between 
model predictions and observed data are needed 
to investigate the star formation histories for late- 
type galaxies. 
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